International  Journal  of  Environment,  Agriculture  and  Biotechnology  (IJEAB) 
http://dx.  doi.  ora/ 10. 221 61/ijeab/2. 5.19 


Vol-2,  Issue-5,  Sep-Oct-  2017 
ISSN:  2456-1878 


Mathematical  Modeling  of  Sun  and  Solar  Drying 
Kinetics  of  F ermented  Cocoa  Beans 

Olabinjo  O.O. v*,  Olajide  J.O.2,  Olalusi  A  .P.1 

'Department  of  Agricultural  Engineering  and  Environmental  Engineering,  Federal  University  of  Technology,  Akure. 
department  of  Food  Science  and  Technology,  LadokeAkintola  University  of  Technology,  Ogbomoso 


Abstract —  In  this  study,  thin  layer  drying  experiments 
were  conducted  to  compute  drying  characteristics  of 
fermented  cocoa  beans  in  open  sun  and  indirect  natural 
convection  solar  dryer.  The  drying  experiments  were 
conducted  at  the  same  time  for  comparison.  Three 
different  thin  layers  drying  of  the  fermented  beans  were 
examined  under  field  conditions  for  Akure,  Nigeria.  The 
drying  process  took  place  only  in  the  falling  rate  period. 
The  drying  curves  obtained  from  the  experimental  data 
were  fitted  to  thirteen  (13)  different  thin  layer 
mathematical  models.  All  the  models  were  compared 
according  to  three  evaluation  parameters.  These  include 
coefficient  of  determination  (R2),  Root  mean  square  error 
(RMSE)  and  Chi-square  (X2).The  results  showed  that 
increasing  drying  air  temperature  resulted  to  shorter 
drying  times.  The  Vermal  et  al.  model  was  found  to  be  the 
most  suitable  for  describing  the  drying  curve  of  the 
convective  indirect  solar  drying  process  of  cocoa  beans 
with  R2  =  0.9562,  X2=0.0069  and  RMSE=0.0067;  while, 
the  Midilli  and  Kucuk  model,  best  described  the  drying 
curve  of  fermented  cocoa  beans  under  open  sun  with  R2  = 
0.9866,  X2 =0.0024  and RMSE=0.0023. 

Keywords —  Thin-layer  drying,  moisture  content, 
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I.  INTRODUCTION 

Drying  is  one  of  the  oldest  methods  of  food  preservation 
(Doymaz,  2007).  Agricultural  and  other  products  have 
been  dried  by  sun  and  wind  in  the  open  air  for  thousands 
of  years.  The  purpose  is  either  to  preserve  them  for  later 
use,  as  in  the  case  with  food;  or  as  an  integral  part  of  the 
production  process  as  with  tobacco  and  cocoa  beans.  It  is 
necessary  that  the  traditional  techniques  be  replaced  with 
industrial  drying  methods.  (Ertekin  and  Yaldiz,  2004). 
Mathematical  modelling  and  simulation  of  drying  curves 
under  different  conditions  is  important  to  obtain  an 
overall  improvement  of  the  quality  of  the  final  product. 
Simulation  models  of  the  drying  process  are  used  for 
developing  new  designs,  improving  existing  drying 
systems,  predicting  the  airflow  over  the  product  and  for 
the  control  of  the  process.  (Aghbashlo,  et  al.,  2008).  Thin 
layer  drying  equations  are  used  to  estimate  drying 


kineticsfor  several  products  and  also  to  generalize  drying 
curves. 

A  critically  important  aspect  of  drying  technology  is 
mathematical  modelling  of  the  drying  process.  Modelling 
of  drying  process  and  kinetics  is  a  tool  for  process  control 
and  necessary  to  choose  suitable  method  of  drying  for 
specific  product.  Modelling  is  also  essential  for  engineers 
to  choose  the  most  suitable  climatic  conditions  in  order  to 
design  appropriate  drying  equipment  for  perishable  crops. 
The  aim  of  this  work  is  to  study  the  drying  process  and 
select  the  most  suitable  model  (in  terms  of  fitting  ability) 
to  describe  the  thin-layer  drying  of  cocoa  beans.  Although 
much  information  has  been  reported  about  modelling  of 
thin  layer  drying  (Togrul&Pehlivan,  2002)  there  is  no 
information  about  modelling  of  thin  layer  drying  of  cocoa 
beans  in  Nigeria. 

II.  MATERIALS  AND  METHODS 

2.1  Drying  experiments 

In  this  study,  fresh  healthy  cocoa  pods  (Amenlonado 
variety)  were  procured  from  Oda  village,  Akure  South 
Local  Government,  Akure  Ondo  State.  The  drying 
experiments  were  carried  out  using  mobile  solar  dryer  in 
the  Department  of  Agricultural  Engineering,  Federal 
University  of  Technology,  Akure.  Plate  1  shows  the 
schematic  diagram  of  the  solar  dryer  used  for  the 
experimental  work  which  consists  of  a  solar  collector  and 
a  drying  chamber.  The  samples  were  weighed  using  a 
digital  balance  with  0.0  lg  sensitivity  every  60  minutes 
throughout  the  drying  process. 

Three  different  thin  samples  of  wet  fermented  cocoa 
beans  were  spread  evenlyinto  the  solar  dryer  and  in  the 
open  sun  drier  for  the  dehydration  test.  The  experiment 
was  replicated  thrice  and  the  mean  value  was  used. 
Thermal  drying  method  was  used  in  the  determination  of 
moisture  content  of  the  samples. lOOg  of  sample  were 
placed  in  oven  at  105+  3°C  and  allowed  to  dry  to  a 
constant  weight  for  24  hours  (Lagha-Benamrouche,  S. 
and  Madani,  K.,  2013).  The  moisture  content  (MC)  was 
calculated  by  expressing  the  weight  loss  upon  drying  a 
fraction  of  the  initial  weight  of  sample  used.  The  moisture 
content  of  the  seeds  was  determined  by  gravimetric 
method  which  determines  the  mass  loss  from  the  sample 
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by  drying  to  constant  weight  (ASABE  STANDARDS, 

1993  and  AOAC,  2000). 

DM (%)  =  *  ioo  (1) 

v  '  w1-w0 

%MCdb  =  100  -  DM%  (2) 

Where  W0  is  weight  of  empty  crucible 


Wtis  weight  of  crucible  plus  sample  before  drying 
144. is  weight  of  crucible  plus  sample  after  drying 
DM  is  dry  matter  and  MCdb  is  the  cocoa  beans  moisture 
content  (g  water/g  dry  base,  d.b). 


Plate.  1:  Mobile  Solar  dryer 


2.2  Mathematical  modelling  of  drying  process 
Many  researchers  have  worked  on  many  thin  layer  models 
in  the  past  and  this  study  evaluate  thirteen  (13)  of  such 
models  as  shown  in  Table  1. 

The  moisture  ratio,  MR  is  given  as  follows: 


MR 


M—Me 

M0~Me 


(3) 


Where  MR  is  the  dimensionless  moisture  ratio  or 
unaccomplished  moisture  content,  M,  Me,  M0  are 
moisture  content  (kg  water/kg,  dry  matter)  respectively. 


The  values  of  Me  are  relatively  small  compared  to  those 
of  M  or  M0  hence  error  involved  in  its  simplification  as 
negligible.  (Aghbashlo,Kianmerhrk&Samini-Akhljahani, 
2008),  hence  moisture  ratio  is  calculated 

(4) 


M  M 

MR  =  —  or  — 

Mr, 


For  drying  model  selection,  drying  data  were  fitted  into 
thirteen  well  known  thin  layer  drying  models  which  are 
given  in  Table  1. 


Table.  1:  Thin  layer  models  used  by  some  researchers  and  used  in  evaluating  the  drying  kinetics  of  Cocoa  beans. 


S/N 

Model  name 

Model  equation 

Refrences 

1 

Newton 

MR=exp(-kt) 

Upadhyayef  al.,  2008 

2 

Page 

MR=exp(-ktAn) 

Saeed  et  al.,  (2006) 

3 

Modified  page 

MR=exp[-(kt)An] 

Ceylan et  al.,  (2007) 

4 

Henderson  and  Pabis 

MR  —  aexp{—kt) 

Kashaninejadand  Tabil(2004) 

5 

Logarithmic 

MR  —  aexp(-kt)  +  c 

Wang  et  al.,  (2007) 

6 

Two-term 

MR=aexp(-ko  t)+bexp(-klt) 

Wang  et  al.,  (2007) 

7 

Two-term  exponential 

MR  —  aexp(— /c0t)bexp(— /qt) 

Tarigan et  al.,  (2007) 

8 

Wang  and  Singh 

MR=l+at+  [bt]  A2 

Wang  and  Singh  (1978) 

9 

Diffusion  approach 

MR=aexp(-kt)+(  1  -a)expi/»i(-kbt) 

Wang  et  al.,  (2007); 

10 

Modified  Henderson  and 

Pabis 

MR=aexp(-kt)+bexp(-gt)+cexp(  -ht) 

Kaya  et  al.,  (2007b) 

11 

Verma  et  al 

MR=aexp(-kt(n))+bt 

Doymaz,  (2005b) 

12 

Midilli  and  Kucuk 

MR  —  aexp^—kt7^  +  bt 

Midillief  al.,  (2002) 

13 

Thomson 

t=a  lnMR+b(In(MR))A2 

Thomson  et  al.,  (1968) 

Moisture  ratio  (MR)  =  dependent  variable,  Drying  constant  constant  (k)  =  independent  variable. 


The  goodness  of  fit  was  determined  using  three 
parameters;  coefficient  of  determination  (R2),  reduced 
chi-square  (x2)  and  the  root  mean  square  error  (RMSE) 


using  equations  (4)  -  (6)  as  in  Sacilik  and  Elicin 
(2008). The  statistical  analyses  were  carried  out  using 
SPSS  13.0  software  and  non-regression  technique. 
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a.  Coefficient  of  determination  (R2) 

R2 

Yfi=1{MRi-MRpre,i)  .Yfl-^MRi-MRexp.i) 
V[Y,f=1(MRi-MRpre,i)\.V\Y,f=1(MRi-MRexp,i)\ 

(5) 

Chi-square  (x2) 

^.2  _  £i=i (MReXp  t~  MRpre.i) 

_  N—n  ’ 

b.  Root  mean  square  error  (RMSE) 

RMSE  =  \^MRVre.i-MRexV.iH/2 

N  J 

Where  MRexp,iis  ith  experimentally  observed  moisture 
ratio,  MRpre,iis  ith  predicted  moisture  ratio,  N  is  the 
number  of  observation,  n  is  the  number  of  model 
constants. 

III.  RESULTS  AND  DISCUSSIONS 


drying  at  any  time  are  represented  in  Figures  1-3.  It  was 
clearly  evident  from  these  curves  that  the  drying  rate  of 
fermented  cocoa  beans  in  the  solar  dryer  was  faster  than 
that  of  the  open  sun  drying.  The  moisture  content  of  the 
fermented  cocoa  beans  reached  6.5%  dry  basis  in  32hours 
of  drying  in  the  solar  dryer,  whereas  the  final  moisture  of 
the  same  product  dried  by  open  sun  drying  was  only 
9.87%  dry  basis  thus  moisture  content  was  not  enough  for 
safe  storage.  When  it  was  dried  under  open  sun  drying, 
the  duration  of  dry  was  about  two  (2)  sunshine  days  to 
bring  it  to  the  same  moisture  level. 

This  can  be  explained  that  the  main  factor  influencing 
drying  rate  was  the  drying  air  temperature.  Compared  to 
open  sun  drying,  solar  dryer  can  generate  higher  air 
temperature  and  affected  the  significant  increasing  of 
evaporation  rate  of  water  and  then  result  in  lower  final 
moisture  content  of  drying  samples.  These  results 
indicated  that  solar  dryer  was  effective  than  open  sun 
drying. 


3.1  Drying  Kinetics  of  Fermented  cocoa  beans 

From  the  experimental  data,  the  moisture  content  (%wb) 

of  fermented  cocoa  beans  for  the  solar  dryer  and  open  sun 
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Fig. 3:  Variation  of  moisture  content  with  drying  time  for  fermented  cocoa  beans  for  2.97g/cm2 


3.3  Mathematical  modelling 

The  moisture  content  data  at  different  experimental 
modes  were  converted  to  the  more  useful  moisture  ratio 
expression,  and  curve  fitting  computations  with  drying 
time  were  performed  with  the  thirteen  (13)  drying  models 
presented  by  previous  workers  (Table  l).The  results  of 
the  statistical  analyses  undertaken  on  these  models  for  the 
natural  convention  solar  drying  and  the  natural  sun  drying 
are  given  in  Table  2  and  3,  respectively.  The  models  were 
evaluated  based  on  Coefficient  of  determination  (R2), 
Chi-square  (x2),  and  Root  Mean  Square  Error  (RMSE). 
All  equations  gave  consistently  high  (R)  values  in  the 
range  of  0.91  -0.98. This  indicates  that  all  equations  could 
satisfactorily  describe  the  solar  drying  rates  of  fermented 
cocoa  beans.  RMSE  ranged  from  0.0067-1.6465,  chi- 
square  ranged  from  0.0069-1.1000  and  for  solar  drying 
while  RMSE  ranges  from  0.0023-1.4672,  and  chi-square 
ranges  from  0.0024-1.2700  for  open  sun  drying.  The 
result  shows  that  for  all  thin  layer  drying  models  and 


conditions  of  solar  drying  ,the  Verma et  al.  (1985)  model 
gave  the  best  fit  with  R2  =  0.9562  ,  x2=  0.0069  ,  and 
RMSE=0.0067  for  solar  drying  .  The  Midilli  model  gave 
the  best  fit  with  R2  =  0.9866,  x2=  0.0024,  and  RMSE  = 
0.0023  for  open  sun  drying.  The  drying  constants  (k)  and 
(1)  and  coefficients  (a)  and  (n)  values  as  well  as  the 
statistical  parameters  R2,  x2,  and  RMSE  are  shown  in 
Tables  2  and  3  for  both  solar  drying  and  open  sun  drying. 
Validation  of  the  Vermal et  al  and  Midilli  and  Kucuk 
models  weremade  by  comparing  the  predicted  moisture 
ratio  with  the  experimented  moisture  ratio  values  from 
all  the  tests.  The  performance  of  the  Vermal  et  al  model 
for  the  thin  solar  drying  and  the  Midilli  and  Kucuk 
model  for  natural  sun  drying  was  illustrated  in  Figures  4 
and  5. The  predicted  data  is  banded  around  the  straight 
line  which  showed  the  suitability  of  the  Vermal  et  al  and 
Midilli  and  Kucuk  models  in  describing  the  drying 
behaviour  of  fermented  cocoa  beans  in  solar  and  open 
drying  respectively. 


Table. 2:  Modelling  the  drying  process  of  fermented  cocoa  beans  using  solar  dryer 


MODELS 

COEFFICIENT 

RA2 

xA2 

RMSE 

NEWTON 

k=0.0186 

0.9311 

0.0134 

0.013 

PAGE 

k=0.0383 

n=0.8255 

0.9351 

0.0108 

0.0105 

MODIFIED  PAGE 

HENDERSON 

k=0.0192 

n=0.8255 

0.9351 

0.0675 

0.0654 

&PABIS 

k=0.0160 

a=0.8748 

0.9506 

0.0085 

0.0082 

LOGARITHMIC 

k=0.0057 

a=  1.5578 

c=-0.7343 

a=0.4523 

0.9662 

0.0088 

0.0082 

TWO  TERM 

ko=0.0162 

kl=0.0160 

b=0.4226 

0.9506 

0.0087 

0.0081 

TWO  TERM  MOD. 

k=0.2539 

a=0.0684 

0.9398 

0.0114 

0.0107 

WANG  &SINGH 

a=-0.0139 

b=5.0149 

0.9234 

1.1 

106465 

APPRO  OF  DIF 

MOD  HENDER 

k=  1.1408 

a=0.1712 

b=0.0133 

0.9562 

0.2088 

0.2021 

&PAB 

k=0.0165 

c=0.3049 

a=0.2746 

g=0.0160 

b=0.2954 

h=0.0158 

0.9506 

0.009 

0.0082 

www.iieab.com  Page  |  2422 


International  Journal  of  Environment,  Agriculture  and  Biotechnology  (IJEAB) 
http://dx.  doi.  ora/ 10. 221 61/ijeab/2. 5.19 


Vol-2,  Issue-5,  Sep-Oct-  2017 
ISSN:  2456-1878 


VERMA  ET  AL 

k=1.1444 

a=0.1705 

g=0.0152 

0.9562 

0.0069 

0.0067 

a=  1.0050, 

MIDILLI  &  KUCUK 

k=0.2021 

n=0.1734 

b=0.0049 

0.9787 

0.4659 

0.4358 

Thompson 

a=-64.7773 

b=-8.0810 

0.9655 

0.8923 

0.9568 

Fig. 4:  Comparison  of  experimental  and  predicted  dimensionless  moisture  ratio  for  solar  drying 


Table. 3:  Modelling  the  drying  process  of  fermented  cocoa  beans  using  open  sun  drying. 


MODELS 

COELLICIENT 

RA2 

XA2 

RMSE 

NEWTON 

k=0.0199 

0.9423 

0.0089 

0.0086 

PAGE 

k=0.0710 

0.9617 

0.0063 

0.0061 

MODIFIED  PAGE 

k=0.0217 

n=0.6905 

0.9617 

0.1376 

0.1339 

HENDERSON 

&PABIS 

k=0.0160 

a=0.8334 

0.9741 

0.0044 

0.0042 

LOGARITHMIC 

k=0.0125 

a=0.9083 

c=-0.0955 

0.977 

0.0038 

0.0037 

ko 

TWO  TERM 

=0.0160 

kl=0.0161 

a=0.4507,b=0.3827 

0.9741 

0.0045 

0.0042 

TWO  TERM  MOD. 

k=0.1656 

a=0.1071 

0.9572 

0.007 

0.0066 

WANG  &SINGH 

a=-0.0144 

b=5.3256 

0.9149 

1.27 

1.1197 

APPRO  OF  DIF 

k=0.8487 

a=0.2269 

b=0.0173 

0.9827 

0.0032 

0.0031 

MOD  HENDER 

&PAB 

k=0.0159 

a=0.2610 

b=0.2814 

0.9741 

0.0048 

0.0042 

c=0.2909 

g=0.0161 

h=0.0161 

VERMA  ET  AL 

k=0.8584 

a=0.0226 

g=0.0147 

0.9827 

0.0135 

0.0132 

MIDILLI  & 

KUCUK 

k=0.1921 

n=0.3235 

a=1.0040,b=-0.0023 

0.9866 

0.0024 

0.0023 

THOMPSON 

a=-59.2342 

b=-4.0169 

0.9827 

0.3246 

1.4672 

www.iieab.com 


Page 


2423 


International  Journal  of  Environment,  Agriculture  and  Biotechnology  (IJEAB) 
http://dx.doi.Org/10.22161/iieab/2.5.19 


Vol-2,  Issue-5,  Sep-Oct-  2017 
ISSN:  2456-1878 


Fig. 5:  Comparison  of  experimental  and  predicted  dimensionless  moisture  ratio  for  open  sun  drying 


IV.  CONCLUSION 

The  solar  dryer  proves  useful  for  local  farmer  as  it  ensures 
high  quality,  good  colour  and  flavour  and  reduces  the 
drying  time.  It  does  not  pollute  the  environment,  requires 
minimal  maintenance  once  it  is  installed  and  with  good 
quality.  In  order  to  explain  the  drying  behaviour  and  the 
mathematical  models  of  fermented  cocoa  beans,  thirteen 
models  were  applied  to  thin  layer  convective  indirect  solar 
dryer  and  open  drying  processes.  The  result  showed  that 
the  Verma et  al  model  was  found  to  be  the  most  suitable 
model  for  describing  the  drying  curve  of  the  convective 
indirect  solar  drying  process  of  cocoa  beans  with  R2 
=0.9562,  jr=0.0068,  MBE=0.0383  and  RMSE=0.0067 
while,  the  Midilli  and  Kucuk  model,  best  described  the 
drying  curve  of  fermented  cocoa  beans  under  open  sun 
with  R2  =  0.9866,  jr=0.0024,  MBE=0.0078  and 
RMSE=0.0023. 
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